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Abstract 

We discuss the strong-coupling expansion in Euclidean field theory. In a formal represen- 
tation for the Schwinger functional, we treat the off-diagonal terms of the Gaussian factor as 
a perturbation about the remaining terms of the functional integral. In this way, we develop 
a perturbative expansion around the ultra-local model, where fields defined at different points 
of Euclidean space are decoupled. We first study the strong-coupling expansion in the (A(/?'^)d 
theory and also quantum electrodynamics. Assuming the ultra-local approximation, we examine 
the singularities of these perturbative expansions, analysing the analytic structure of the zero- 
dimensional generating functions in the complex coupling constants plane. Second, we discuss 
the ultra-local generating functional in two idealized field theory models defined by the following 
interaction Lagrangians: Cxx{gi, g2] (p) = giip^{x) + g2ip~''^{x), and the sinh-Gordon model, i.e., 
Cxxx{g3, g^] tp) = g^ {cosh.{g4^ip{x)) — 1). To control the divergences of the strong-coupling pertur- 



bative expansion two different steps are used tfirougfiout tlie paper. First, we introduce a iattice 
structure to give meaning to tfie ultra-local generating functional. Using an analytic regulariza- 
tion procedure we discuss briefly how it is possible to obtain a renormalized Schwinger functional 
associated with these scalar models, going beyond the ultra-local approximation without using a 
lattice regularization procedure. Using the strong-coupling perturbative expansion we show how 
it is possible to compute the renormalized vacuum energy of a self-interacting scalar field, going 
beyond the one-loop level. 
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1 Introduction 



The purpose of this paper is twofold. The first is to study the structure of the singularities of the 
perturbative series in different perturbatively renormalizable models in field theory. Our method is 
based in the fact that these non-perturbative results can not be obtained using the weak coupling 
perturbative expansion, and a different perturbative expansion is mandatory. Thus, in the context 
of the strong-coupling perturbative expansion, we investigate the analytic structure of the ultra- 
local generating functionals for the {X<fi'^)d model and quantum electrodynamics in the complex 
coupling constant plane. The second is to discuss how the strong- coupling perturbative expansion 
can be used as an alternative method to compute the renormalized free energy or the vacuum 
energy associated with self-interacting fields, going beyond the one-loop level. To implement 
such method we have to show that it is possible to obtain renormalizable Schwinger functionals 
associated with scalar models going beyond the ultra-local approximation. This is performed using 
an analytic regularization procedure. This work is a natural extension of the program developed 
by Klauder [1], Rivers [2] and others who have been studying the strong-coupling expansion and 
the ultra-local generating functional in different scalar infrared-free models in field theory. 

The perturbative renormalization approach in quantum field theory is an algorithm where, 
starting from the Feynman diagrammatic representation of the perturbative series, two different 
steps are usually performed. In the first step, we control all the ultraviolet divergences of the 
theory, using a procedure to obtain well-defined expressions for each Feynman diagram. In the 
second step we have to implement a renormalization prescription where the divergent part of 
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each Feynman diagram is canceled out by a suitable counterterm. For complete reviews of this 
program see for example Ref . [3] or Ref . [4] . Concerning the first step, there are different ways to 
transform the Feynman diagrams in well-defined finite quantities. The simplest way is to modify 
the field theory at short distances by introducing a sharp cut-off in momentum space, or a more 
elaborate version, as the Pauli-Villars regularization [5] . In this second case we simply modify the 
propagator for large momentum. We can also employ a lattice regularization method in which 
we replace the continuum Euclidean space by a hypercubic lattice, with lattice spacing a. It is 
clear that the introduction of a lattice provides a cut-off in momentum space of the order of the 
inverse of the lattice spacing a. Finally, a more convenient regularization procedure is dimensional 
regularization [6] [7] [8] [9], which is particulary well suited to deal with abelian and non-abelian 
gauge theories. 

In the second step, to implement the renormalization procedure, using dimensional regular- 
ization for example, we use the fact that the ultraviolet divergences of Feynman diagrams appear 
as poles of some function defined in a complex plane. Then the perturbative renormalization is 
performed by the cancelation of the principal part of the Laurent series of the analytic regularized 
expressions. This cancelation is done introducing counterterms in the theory. There are different 
ways to disregard the divergent part of each Feynman diagram. For example, we can use the 
minimal subtraction scheme (MS) where the counterterms just cancel the principal part of the 
Laurent series of the analytic regularized expressions, or any different renormalization scheme. 
The arbitrariness of the method employed must be cured by the renormalization group equations 
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[10] [11] [12]. 

In this framework, field tfieory models are classified as perturbatively super-renormalizable, 
renormalizable and non-renormalizable. In the standard weak-coupling perturbative expansion, 
the fundamental difference between a perturbatively renormalizable and a non-renormalizable 
model in field theory is given by the following property. In a non-renormalizable model, the 
usual renormalization procedure used to remove the infinities that arise in the usual perturbative 
expansion introduces infinitely many new empirical parameters in the theory. This comes from 
the fact that we need to specify the finite part of infinitely many counterterms. Consequently, 
the predictibility or the physical consistency of non-renormalizable models is missing. Thus, 
renormalizability of field theory models provides a valuable constraint on new theories. See, for 
instance, the discussion in Ref. [13]. Nevertheless, there are some non-renormalizable models 
where it is possible to construct a physically sensible version of the theory. A well-known example 
is the Gross-Neveu model, which is not renormalizable in the usual sense. However, this model 
is renormalizable in the expansion, for < 4 [14] [15]. We would like to mention that there 
is a alternative approach to deal with non-renormalizable model in field theory, using the idea of 
effective field theories [16] [17], but in this paper we will not discuss these issues. 

In the case of perturbatively renormalizable models, although the renormalization procedure 
can be implemented in a mathematicaly consistent way, it is still not clear how the renormalized 
perturbative series can be summed up in different models [18] [19] [20]. In the literature, there 
are many results showing that the series that we obtain in different perturbative renormalizable 
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theories in o? = 4 do not converge for any value of the couphng constants of the interacting theories. 
Well-known theories with such problems are scalar models with a (A<^^) self-interaction and also 
quantum electrodynamics. If one tries to perform a partial resummation of the perturbative 
series in both theories, Landau poles appear [21] [22] [23]. In a four-dimensional spacetime, to 
circumvent the problem of non-hermitian Hamiltonians in the infinite cut-off limit, the (A(/?^)4 
model, the 0(7V)4 model, and also quantum electrodynamics must be made trivial. 

A new step in the development of quantum field theory was given by the construction of 
non-abelian gauge field theories and the discovery of asymptotic freedom [24] [25] [26] after the 
construction of the renormalization group equations. From the renormalization group equations 
an important classification of different field theory models arises. The models are either asymptot- 
ically free or IR (infrared) stable. In the renormalization group approach, the triviality of {Xlp'^)4 
model, the 0{N)4 model and also quantum electrodynamics in a four-dimensional spacetime is a 
refiection of the absence of a non-trivial ultraviolet stable fixed point in the Callan-Symanzik /3- 
function. In the infrared-free theories, for d — 4, the problem of the singularities of the connected 
three-point and four-point Green's functions for quantum electrodynamics and (A</?'^)4, respec- 
tively, is related to the following fact: in the framework of the weak-coupling expansion, the high 
frequency fluctuations are more strongly coupled than the lower frequency ones. It is well known 
that we found a completely distinct behavior in a theory where the high frequency fluctuations are 
more weakly coupled than the lower frequencies ones. In this situation, at least the zero charge 
problem does not appear. 
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The unification of statistical meclianics witli some models in quantum field theory was achieved 
in progressive steps. First Schwinger introduced the idea of Euclidean fields, where the classical 
action must be continued to Euclidean time [27]. Then, Symanzik constructed the Euclidean 
functional integral where the vacuum persistence functional, defined in Minkowski spacetime, 
becomes a statistical mechanics average of classical fields weighted by a Boltzmann probability 
[28] [29]. At the same time, a deeper insight into our understanding of the renormalization 
procedure in different models in field theory was given by the study of the critical phenomena and 
the Wilson version of the renormalization group equations [30]. Further, Osterwalder and Schrader 
proved that for scalar theories, the Euclidean Green's functions or the Schwinger functions, which 
are the moments of the Boltzmann measure, are equivalent to the Minkowski Green's functions 
[31]. A new step in our understanding of the limitation of the perturbative approach in quantum 
field theory was achieved by Aizenman [32] and Frohhch [33]. These authors proved that the 
{X(f'^)d model, with the use of a lattice regularization with nearest neighborhood realization of the 
Laplacian, leads to a trivial theory in the continuum limit for (i > 5. For d — 4, with additional 
assumptions, it is also possible to obtain the triviality of the model. For an interesting review of 
the triviality problem in quantum field theory, see Ref. [34]. 

It is important to point out that some authors claim that the triviality of {Xip^)d for > 4 is 
an odd result, since for d — 4 the renormalized perturbative series is non-trivial and for (i > 5 the 
theory is perturbatively non-renormalizable. For example, Klauder has argued that the triviality 
of {X(p'^)d for (i > 4 is still an open problem [35]. Making use of the correspondence principle, this 
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author has emphasized that the quantization of a non-trivial classical theory can not be a non- 
interacting quantum theory. Furthermore, he claims that an alternative regularization procedure 
can give a non-trivial theory in the infinite cut-off limit. 

Some results going in this direction have been obtained by Gallavotti and Rivasseau [36]. 
These authors discuss the scalar {Xy:>'^)d theory with more general regularized theories where the 
realization of the Laplacian is not restricted to the nearest neighbours and also the presence of 
antiferromagnetic couplings. They suggested that the ultraviolet limit of such lattice regularized 
field theories is not a Gaussian field theory model, which would open the possibility to construct 
scalar models with a non-trivial ultraviolet limit. We shall emphasize that the ultraviolet behavior 
of the {X(fi'^)d model for d > 4 is a strong-coupling problem, since, in the weak coupling perturbative 
expansion, the high-order terms of the perturbative series are dominant in the large cut-off limit, 
as has been discussed by many authors. See, for instance, the discussion in Ref. [37]. Another 
situation in field theory where the weak-coupling perturbative expansion is not appropriate is in 
the large distance behavior of quantum chromodynamics. 

In the case of a strong-coupling regime of a theory, we can try to obtain rigorous results by 
the use of constructive field theory, perform a partial resummation of the Feynman diagrams 
improving the Feynman-Dyson perturbative series, or perform a different perturbative expansion 
by using the following approaches. The first is to introduce auxiliary fields in order to disconect the 
interaction part from the free part of the Lagrangian density and then to perform a perturbative 
expansion of the Schwinger functional in inverse powers of the coupling constant [38] . For example. 
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the expansion, where N is the number of the components of the field in some isotopic space 
(the dimension of the order parameter), is a reahzation of this approach [39]. Of course we are 
still using the standard perturbative scheme, performing a perturbative expansion with respect 
to the anharmonic terms of the theory. The basic idea of the second approach is to construct a 
formal representation for the generating functional of complete Schwinger functions of the theory, 
treating the off-diagonal terms of the Gaussian factor as a perturbation about the remaining terms 
in the functional integral. This approach has been called the strong-coupling expansion [40] [41] 
[42] [43] [44] [45]. 

The purpose of this paper is to discuss some of the problems and virtues of the strong-coupling 
expansion in different scalar models and also quantum electrodynamics. We first study the singu- 
larities of the generating functional of complete Schwinger functions for the {XLp^)d model in the 
complex coupling constant plane by examining the analytic structure of the zero-dimensional gen- 
erating function. Then, we repeat the analysis for the case of quantum electrodynamics. Second 
we present two idealized interacting field theory models where the weak-coupling expansion cannot 
be used. In these cases the strong-coupling expansion is more adequate to investigate the proper- 
ties of the models. It is important to stress that the analytic structure of the Schwinger functions 
(in the weak-coupling perturbative expansion framework) and also of the Schwinger functional 
can be easily obtained using the ultra-local approximation derived in the strong-coupling pertur- 
bative expansion. As we will see, these non-perturbative results do not change if we go beyond 
the ultra-local approximation. We should say at this point that Bender et al [46] also consider 
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the zero-dimensional field theory to obtain non-perturbative results in field theory. Third, using 
an analytic regularization procedure, we discuss briefly how it is possible to obtain a renormal- 
ized Schwinger functional going beyond the ultra-local approximation. Finally, we sketch how 
the strong-coupling perturbative expansion can be used as an alternative method to compute the 
renormalized free energy or the vacuum energy associated with self-interacting fields, going beyond 
the one-loop level. 

There are three points that we would like to briefly discuss. First, it is interesting to known 
whether an ultra-local approximation has been used in the literature in other contexts. In the 
Laudau theory of continuous phase transition, which reproduces the mean-fleld exponents, we have 
a simplifled version of the ultra-local approximation, since in the partition function we drop the 
gradient term and the sum will be dominated by its largest term. Further, Landau and Ginzburg 
modifled the the original Landau theory by introducing the gradient term into the energy density 
that discourages rapid fluctuations in the order parameter [47]. Before continuing we would 
like to point out that the strong-coupling perturbative expansion is quite similar to the high- 
temperature series expansion in statistical mechanics. A similar idea is the Mayer expansion, a 
method for carrying out the cluster expansion, introduced into quantum fleld theory by Symanzik 
[48] . Finally, in the study of critical phenomena using lattice simulation, there is an analog of the 
strong-coupling perturbative expansion: the hopping parameter expansion, where the perturbative 
expansion starts from the disordered lattice system [49] . For a complete review of the linked cluster 
expansion, see for example Ref. [50]. 
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Second, it is interesting to point out that the study of the ultra-local field theory models can 
also bring some insights into non-renormalizable field theory models. As a simple example, let us 
suppose a massive abelian vector field Wi^{x) coupled with fermions. Using a Fourier representation 
for the Euclidean two-point Schwinger function associated with the massive vector field 
the Fourier transforms Gi^i,{k) are given by 



The Fourier transform G^y{k) of the two-point Schwinger function does not vanish in the ultraviolet 
limit k ^ oo and its large Euclidean momentum behavior is roughly m~^. Consequently the 
behavior in this limit is similar to the behavior for the two-point Schwinger function in the ultra- 
local {\ip'^)ii model. 

Finally, we would like to call the attention of the reader that the study of the analytic structure 
of theories in the complex coupling constant plane has been used by many authors in quantum 
field theory. It is well known that the behavior of the standard perturbative series in powers of 
the coupling constant at large order is related to the analytic structure of the partition function in 
a neighborhhood of the origin in the complex coupling constant plane. For example. Bender and 
Wu [51] studied the anharmonic oscillator and pointed out that there is a relation between the 
n^^ Rayleigh-Schrodinger coefficients and the lifetime of the unstable states of a negative couphng 
constant anharmonic oscillator. 

The organization of the paper is as follows: In section II we discuss the standard weak-coupling 
expansion for the {\ip'^)d model. In section III we discuss the strong-coupling expansion for the 
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(A(/?^)d model. In section IV we perform the study of the analytic structure of the zero- dimensional 
(A<^^) model in the complex coupling constant plane. In section V we perform the study of 
the analytic structure of the zero-dimensional quantum electrodynamics in the complex coupling 
constant plane. In section VI we present idealized models where the strong-coupling perturbative 
expansion must be used. In section VII we sketch how it is possible to go beyond the ultra- 
local approximation in general scalar models. In section VIII we show how it is possible to 
compute the renormalized vacuum energy associated with a self-interacting scalar field in the 
presence of macroscopic structures going beyond the one-loop level. Finally, section IX contains 
our conclusions. To simplify the calculations we assume the units to be such that h = c = 1, and 
also all the physical quantities are dimensionless. Consequently it is convenient to introduce an 
arbitrary parameter // with dimension of mass to define all the dimensionless physical quantities. 
Thus, in the paper we are using dimensionless cartesian coordinates (x" = fj,x"^), where x"^ are 
the usual coordinates with dimension of length. 

2 Weak coupling perturbative expansion for the scalar 
{X(f^)d model 

Let us consider a neutral scalar field with a (A(/?^) self-interaction, defined in a d-dimensional 
Minkowski spacetime. The vacuum persistence functional is the generating functional of all vac- 
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uum expectation value of time-ordered products of the theory. The Euchdean field theory can be 
obtained by analytic continuation to imaginary time supported by the positive energy condition 
for the relativistic field theory. In the Euclidean field theory, we have the Euclidean counterpart 
for the vacuum persistence functional, that is, the generating functional of complete Schwinger 
functions. The {Xf'^)d Euclidean theory is defined by these Euclidean Green's functions. The 
Euclidean generating functional Z{h) is formally defined by the following functional integral: 

Z{h) = J [dip] exp (^-So -Si + J d'^x h{x)ip{x)^ , (1) 

where the action that usualy describes a free scalar field is 

So{^) = I d'x (^^{d^r + \ml ^\x)) , (2) 

and the interacting part, defined by the non-Gaussian contribution, is 

S:{^)= jd'x^j^^^'ix). (3) 

In Eq.(l), [dip] is a translational invariant measure, formally given by [df] = JJ^ dp{x). The terms 
go and mg are respectivelly the bare coupling constant and mass squared of the model. Finally, 
h{x) is a smooth function that we introduce to generate the Schwinger functions of the theory by 
functional derivatives. Note that we are using the same notation for functionals and functions, for 
example Z{h) instead the usual notation Z[h]. 

In the weak-coupling perturbative expansion, which is the conventional procedure, we perform 
a formal perturbative expansion with respect to the non-Gaussian terms of the action. As a 
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consequence of this formal expansion, all the n-point unrenormalized Schwinger functions are 
expressed in a powers series of the bare coupling constant gg. Let us summarize how to perform the 
weak-coupling perturbative expansion in the {X(p'^)d theory, and briefly discuss also the divergent 
behavior of the perturbative series. The Gaussian functional integral Zo{h) associated with the 
Euclidean generating functional Z{h) is 

Zo{h) = N j [di^] exp K + h^^ . (4) 

We are using the compact notation of Zinn- Justin [52] and each term in Eq. (4) is given by 

<pKip^ j d'^x j d%(p{x)K{mo; x,y)ip{y), (5) 

and 

hcp = j d'^xip{x)h{x). (6) 
As usual A/" is a normalization factor and the symmetric kernel K{mo] x, y) is deflned by 

K{mo]x,y) = {-A + ml)6\x-y), (7) 

where A denotes the Laplacian in R^. As usual, the normalization factor is defined using the 
condition Zo{h)\h=o — 1. Therefore M — [det(— A -\-m^)Y but, in the following, we are absorbing 
this normalization factor in the functional measure. It is convenient to introduce the inverse kernel, 
i.e. the free two-point Schwinger function GQ{mQ]x — y) which satisfies the identity 

j d'^z Go{mo; X - z)K{mo; z - y) =S'^{x-y). (8) 
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Since Eq.(4) is a Gaussian functional integral, simple manipulations, performing only Gaussian 
integrals, give 

J [dip] exp (^-So + J d'^x h{x)ip{x)^ = exp (^^ j d'x J d'^y h{x) Go{mo; x - y)h{y)^ . (9) 

Therefore, we have an expression for Zo(/i) in terms of the inverse kernel Go{m.o;x — y), i.e., in 
terms of the free two-point Schwinger function. 

This construction is fundamental to perform the Feynman-Dyson weak-coupling perturbative 
expansion with the Feynman diagramatic representation of the perturbative series. Using Eqs.(l), 
(2), and (9), we are able to write the generating functional of all bare Schwinger functions Z{h) 



as 



Z{h) = exp ^- J d^x cxp Q J d^x j d^y h{x) Go{mo; x - y)h{y)) , (10) 

where Cx is defined by the non-Gaussian contribution to the action. Consequently, in the con- 
ventional perturbative expansion, the generating functional of complete Schwinger functions is 
formally given by the following infinite series 

Z{h) = E ( / d^x J^)) exp Q j d'x J d'^y h{x) G^{m^- x - y)h{y)^ . (11) 

To generate all the n-point Schwinger functions we have only to perform a suitable number of 
functional differentiations in Z{h) with respect to the source h{x) and set the source to zero in 
the end. The bare n-point Schwinger functions are defined by 

5 t 



Gn(xi,X2,..,Xn) ^ Z ^(/l = 0) 



5h{xi) 6h{xn) 



m 



\h=0- 



(12) 
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This general method can be used to derive the weak-couphng perturbative expansion in different 
theories. Observe that it is possible to generalize this formalism including the product of composite 
sources [53] [54] [55], but in this paper we limit ourselves to models without composite operators. 
To generate only the connected diagrams G!if>{xi, ..,Xn), let us consider the generating functional 
of the connected Schwinger functions (the free energy functional), defined as F{h) = In Z{h). The 
functional Taylor expansion of the generating functional of connected Schwinger functions is 

oo 1 „ n n 

nh) = E - / n d'x, n gi^h^i, -,^n). (13) 

n=l •' k=l k=l 

The functional Taylor expansion of the generating functional of all Schwinger functions Z{h) in 
powers of the coupling constant go is 

Z(h)=M(l + f:g'oAk{h)y (14) 

where A^^h) are perturbative coefficients. After a regularization and renormalization procedure it 
is possible to show that any physically measurable quantity f{g) can be expanded in power series 
defined by 

00 

f{9) = j:fk9^ (15) 

k=0 

where g is the renormalized coupling constant and fk are perturbative coefficients. 

It is well known that in general the series that we obtain from perturbatively renormalizable 
theories are divergent. Consequently the perturbative renormalization method is not enough to 
obtain well defined physical quantities. For example, for P{(p)2 the renormalized perturbative 
series for any connected Schwinger function that can be obtained by a Wick ordering is divergent 
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[57]. For the {X(p'^)3 model a similar divergent behavior was proved by de Calan and Rivasseau 
[58] . Although in general the series that we obtain from perturbatively renormalizable theories are 
divergent, there is an impressive agreement of the theoretical results with the experiments, when 
someone use the first terms of the weak-coupling perturbative expansion to extract predictable 
results. Therefore, these divergent series shall be an asymptotic expansion of the solutions of the 
theories. In other words, in a specific theory even though the renormalized perturbative series 
diverges, a finite number of terms of the series is still a good approximation of the functions in 
question. Now we need a tool for obtaining the functions of the theory from the divergent series. 
The Borel resummation is this tool that allows us to obtain the solutions of the theory from these 
divergent series [59] . Let us discuss with more details the asymptotic expansion of a function and 
Borel summability. Consider a function f{z) defined in the complex plane for large z. The formal 
series J2'^=o Q^n^~", which need not converge for any value of z, is called the asymptotic expansion 
or asymptotic representation of the function f{z) if defining 

N 

i) Sn{z) = E «n-s~" (16) 

n=l 

It) RN{z) = z^\f{z)-SN{z)\, (17) 

we have lim^zi^^ Rj^{z) = 0, for every fixed A^. There is a similar definition of the asymptotic 
expansion of a function near zero, involving series of the kind J2'^=o <^nZ"'- Again, this series is the 
asymptotic representation of f{z), i.e., f{z) ~ Z^^o '^n^"' if for a small z we still have (i) and in 
(ii) we have Rn{z) — z~^\f{z) — Sn{z)\ and for every fixed N we have lim\z\-^oRN{z) — 0. Prom 
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the above definition there are two main questions. The first is the question whether a function 

under consideration possesses an asymptotic expansion, which we call the expansion problem. 

There is also the question of how the function is to be found, which is represented by a given 

asymptotic expansion, that we call the summation problem. Note that any function can have 

only one asymptotic expansion, or we can show that the function in question has no asymptotic 

expansion. Suppose f{x) — exp{—x), x > 0. It is clear that all the coefficients of the asymptotic 

expansion are zero since a;^exp(— x) — > for every A; > when x ^ oo. This simple result shows 

that different functions may have the same asymptotic expansion. For example if h{z) has an 

asymptotic representation it is clear that h{z) + exp(— 2;) or h{z) + aexp{—bz) for 6 > have 

the same asymptotic representation. Now, consider a function f{z) which has the asymptotic 

expansion in a region of the complex plane, defined by a divergent series. Thus we have 

00 

f{z)r^^hz'. (18) 

fe=0 

The Borel transform of f{z), called Bf{z), is defined as 

00 -1 

Bfiz) = T.T:^h^'- (19) 

The key point is that the Borel transform may converge even if the series is divergent. The Borel 
resummation of the series is obtained applying the inverse Borel transform on Bf{z), given by 

roo 

f{z) = / exp(-t) Bf{zt) dt. (20) 

J 

This construct is an indispensable tool to recover a function from its asymptotic expansion in 
quantum field theory. A pedagogical discussion of the application of the Borel transform in 
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perturbation theory can be found in Ref. [60]. It is not difficult to repeat this construction for 
the n-point Schwinger function of any renormahzable theory. If the perturbative series of the 2n- 
point renormalized Schwinger functions does not converge, as these series must be an asymptotic 
expansion for the solutions of our theory, the Borel resummation method can be used to recover the 
solutions. In the standard perturbative expansion we express the 2n-point renormalized Schwinger 
functions as the following power series in the renormalized coupling constant g: 

oo 

G2n(9; Xi,X2,..,X2n) ^ ^ GS (^1^ ^2, X2n) ■ (21) 

fe=0 

Let us define the Borel transform of the n-point Schwinger function by 

oo 

G'2ji(t; Xi,X2, ..,X2n) = X] TJ G2n i^l, X2, ■ ., X2n) , (22) 

fe=0 

and it is clear that from the inverse Borel transform we have 

I roo j 

G2n{g; Xi,X2, ..,X2n) = - (ir exp( )G2n{r; Xi,X2, ..,X2n)- (23) 

g Jo g 

The series which defines G2n{T; xi,X2, ..,X2n) in Eq. (22) has much better convergence proper- 
ties than the original series of G2n{g] Xi,X2, ■■,X2n) in Eq.(21). Eckmann et al. [61] obtained the 
Schwinger functions of the {Xf*)2 model from the divergent perturbative series using Borel re- 
summation method. The extension of this result to (A<^^)3 model was obtained by Magnen and 
Seneor [62]. 

Note that there are some situations where the Borel resummation method can not be im- 
plemented. This happens, for example, when the Borel transform has singularities in the real 
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line. These singularities are related to non-perturbative effects which are not apparent in the 
weak-coupling perturbative expansion [63]. In the absence of singularities in the positive axis of 
the Borel transform, the Borel resummation method is a powerful way to extract results from a 
divergent series. See also Ref . [64] and Ref . [65] . 

3 The strong-coupling perturbative expansion for the scalar 
{X(p^)d model 



Many phenomena in quantum field theory cannot be described in the framework of the weak- 
couphng perturbative expansion. For example, in the case of the strong-coupling regime of a 
theory, the perturbative expansion in powers of the coupling constant is unrealiable. From now 
on, we shall study the strong-coupling perturbative expansion in different models in field theory. As 
we discussed this perturbative expansion may be used in the strong-coupling regime of a model, as 
for example the ultraviolet limit of a non-asymptoticaly free model or as an alternative expansion 
in models with non-polynomial Lagrangian interaction. In the first situation, or in the case of 
non-renormalizable theories, it is imperative to investigate an alternative perturbation expansion. 
Before continue, we would like to call the attention of the reader that a alternative perturbative 
program for dealing with non-renormalizable theories, has been developed by Klauder and others 
[67] [68] [69] [70]. Klauder proposed a non-canonical formulation for the quantization of (A(/?^)4 
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model {p > 4) using a non-translational invariant functional measure. To support this approach he 
observes that there are many situations where an infinitesimal perturbation causes a discontinuous 
change in the eigenfunctions and eigenvalues associated with a Hamiltonian system [71] [72]. For 
example in the (A<^)4 model, if p < 4 the theory is renormalizable and field configurations which 
have a finite free action also give a finite contribution for the interaction term. For p > 4 this 
does not happen: the free field situation can not be obtained when A — > 0"^ and this limit is the 
pseudo-free solution. In this paper we decided to follow a more conventional treatment. We use 
an alternative perturbative expansion, but assuming that the measure in the functional integral 
is translational invariant instead of using a scale covariant measure. 

The lesson that we have from these discussions is that if someone decides to perform a perturba- 
tive expansion of a strongly-coupled theory, then a resummation of the weak-coupling perturbative 
series to obtain non-perturbative results is necessary. An alternative procedure is not to use the 
conventional perturbation theory around the Gaussian-free theory. Consequently, we now turn to 
the alternative expansion that has been called in the literature the strong-coupling perturbative 
expansion. The basic idea of this approach is to treat off-diagonal terms of the Gaussian factor 
as a perturbation about the remaining diagonal terms in the integral. Although we are studying 
only scalar models, the extension to higher spin fields is straightforward. See for example the 
discussion given in Ref. [73]. 

Let us suppose a compact Euclidean space with or withouth a boundary. An equivalent 
possibility is to work in an unbounded Euclidean space but assume that the functional integral is 
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taken over field configurations that vanish at large Euclidean distances. Let us suppose that there 
exists an eliptic, semi-positive, and self-adjoint differential operator D acting on scalar functions 
on the Euclidean space. The usual examples are D = (—A) and D = (—A + mg). The kernel 
K{mQ] x — y) is given by K{mQ] x — y) — D 5'^{x — y). Let us study first the self-interacting (At^^)^ 
model. Thus we have: 

C^{go;v)^%v\x). (24) 

Treating the off-diagonal terms of the Gaussian factor as a perturbation about the remaining terms 
in the integral, we get a formal expression for the generating functional of complete Schwinger 
functions Z{h): 

=exp(-i/ A /<i^,^^(„„;.-,)^) Q„(„), (26) 

where the ultra-local generating functional Qo{h) is defined by 

Qo{h) =N j [di^] exp (I d^x ( -^-^_'^\x) + h{x)i^{x))^ . (26) 

The factor A'' is a normalization that can be found using that (5o(/i)U=o = 1- Observe that all the 
non-derivative terms in the original action appear in the functional integral that defines Qo{h). 
As we discussed in section II, in the {Xf'^)d model, the kernel is defined by Eq.(7): K{mo; x — y) = 
(—A -I- mg) 5'^{x — y). At this point it is convenient to consider h{x) to be complex. Consequently 
h{x) = Re(/i) -|- ilm{h) (we are concerned with the case Re(/i) = 0). It should be noted that 
although the functional integral Qo{h) is not a product of Gaussian integrals, it can be viewed 
formally as an infinite product of ordinary integrals, one for each point of the d-dimensional 
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Euclidean space. The fundamental problem of the strong-coupling expansion is how to construct 
non-Gaussian measures to define the Schwinger functional. It is important to point out that the 
solution of this problem would allow us to deal with non-renormalizable models in the weak- 
coupling expansion framework. 

The expansion of the exponential term on Eq.(25) gives a formal expansion of Z{h) as a 
perturbative series in the following form: 

oo 

Zih)^^Z^\h), (27) 

i=0 

where the two first terms of the perturbative series are respectively the ultra-local generating 
functional Qo{h) and Z^^\h) defined by 

= (4/ A /<i^,^A-(„„;.-,)^) (28) 

The main difference from the standard perturbative expansion is that we have an expansion of the 
generating functional of complete Schwinger functions in inverse powers of the coupling constant. 
We are developing our perturbative expansion around the static ultra-local functional Qoih) [66] 
[74] [75] [76] . Fields defined in different points of the Euclidean space are decoupled in the ultra- 
local approximation since the gradient term is dropped. 

As we stressed, although the strong-coupling expansion is very inconvenient for pratical cal- 
culations in the continuum Euclidean space, it is very natural in the lattice. The technical 
problems that we have to deal with in the continuum Euclidean space are the following: first we 
have to define non-Gaussian functional measures; second we have to regularize and renormalize 
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the Schwinger functions obtained from the generating functional, going beyond the ultra-local ap- 
proximation. We would like to stress that we will not use a lattice structure of the Euclidean space 
as a regulator to implement the renormalization program. Instead, we use the lattice structure 
only to define what we mean by the ultra-local generating functional Qo{h). 

In this paper we are not interested in regularizing the complete series of the strong-coupling 
perturbative expansion. We are interested first, in the analytic structure of the Schwinger func- 
tional in the complex coupling constant plane. We will show that the zero- dimensional generating 
function of the X(fi'^ model has a branch point singularity. This kind of singularity in the complex 
coupling constant plane is also present in the ultra- local quantum electrodynamics [77]. Our sec- 
ond interest is to study scalar models where the use of the strong-coupling expansion is mandatory. 
Third, we briefly discuss how it is possible to use an analytic rcgularization procedure, in order to 
obtain a renormalized Schwinger functional, going beyond the ultra-local approximation. Finally, 
we sketch how the strong-coupling perturbative expansion can be used as an alternative method 
to compute the renormalized free energy or the vacuum energy associated with self-interacting 
fields, going beyond the one-loop level. 

Using the fact that the functional integral which defines Z{h) is invariant with respect to the 
choice of the quadratic part, let us consider a slightly modification of the strong-coupling expansion. 
We split the quadratic part in the functional integral proportional to the mass squared in two 
parts: the off-diagonal terms of the Gaussian factor and the ultra-local generating functional. The 
new generating functional of the complete Schwinger functions will be defined by the following 
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functional integral 

ZW=exp(-l/d''x/<i^!,^A-(m„..;x-!,)^) Q„(a; k). (29) 

where Qoio", h), the new ultra- local functional integral, is given by 

Qo{<j;h) ^Afjidif] exp (^J d'^x {—amlip^x) - ^^^\x) + h{x)ip{x))j . (30) 

and the new kernel K{mo, a;x — y) is defined by 

K(mo, a;x-y)^ (-A + (1 - a)ml ) S'^ix - y), (31) 

where cr is a complex parameter defined in the region < Re{a) < 1. The choice of a suitable a 
will simplify our calculations in some situations. 

In the next section we study the analytic structure of the ultra- local {X(p'^)d theory. The 
divergences that appear in the formal representation of the Schwinger functional Z{h) are of two 
kinds. The first kind is related to the infinite volume and continuum hypotesis of the Euclidean 
space and this divergence can be controled by the introduction of a box and a regulator function. 
The second kind of divergences is related to the functional form of the non-Gaussian part of the 
action and appears as a divergent perturbative series. We are first concerned with this second 
kind of divergences. 

The study of the ultra-local model in different field theories will clarify the structure of the 
singularities for the perturbative series, and also the structure of the singularities for the n-point 
Schwinger functions in the complex coupling constant plane. We first investigate the analytic 
structure of the {X<^^)d ultra-local model, in the complex coupling constant plane. 
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4 The analytic structure of the ultra-local {X(p'^)d model. 



The aim of this section is to analyze the analytic structure of the zero-dimensional X(p^ model 
in the complex coupling constant go plane. As we discussed before, the first term of the strong- 
coupling expansion of Z(h) is exactly the ultra-local model [66] [74] [75] [78] [79], also called 
the static independent value model. Since we can interpret the ultra-local model as an infinite 
product of ordinary integrals, let us introduce a Euclidean lattice and analyse the generating 
function defined in each point of the Euclidean lattice given by 



where for simplicity we are assuming a — 1. The generating function in the absence of external 
source is defined by 2; (mo, go] h)\h=o = 2:0(^0, go)- Consequently, our aim is to analyse the following 
integral with a quartic probability distribution in which the zero-dimensional partition function 

^o("io,5'o) is given by 



Note that this integral is well defined for Re go > 0. As the exponential power series is 
convergent everywhere, we may write 




(32) 




(33) 




(34) 
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The above series is not uniformly convergent. Therefore we can not interchange the integration 
and the summation. Nevertheless let us perform this interchange, integrate term by term, and get 
a formal series z^^^mo, go). This formal series is the asymptotic expansion for zo{mo,go). Thus 
we have zo{mo,go) ~ z^^\mo,go) and choosing ml = 1 it is not difficult to show that 

oo 

z^^\rnQ, go)\ml=i = Yl i-9o)''ck, (35) 

fe=0 

where the coefficients Cfc are given by = ^^^jf^- The formal series z^^\nio,go) is divergent for 
any non-nule value of go. The asymptotic expansion for the zero- dimensional partition function 
given by z^^\m.o, go) has the contribution from the vacuum diagrams [80] [81], and each coefficient 
Ck is given by the sum of symmetry factors over all diagrams of order k. 

It is not necessary to go so far. In this particular model we will take a short cut. The integral 
given by Eq.(33) can be solved exactly for Re go>0, yielding 

,„).(_L)}™5>fg_|^M). (36) 

where 'i/{a,c; z) is the confluent hypergeometric function of second kind [82] [83], and we are using 
the principal branch of this function. Since we are interested in studying the analytical structure 
of zo{mo,go) in the coupling constant complex plane at ^-o = 0, we must investigate the analytic 
structure of the ^{a,c; z) at z — oo. We are following part of the discussion developed in Ref. [84]. 
The confluent hypergeometric function of second kind ^'(a, c ; z) is a many valued analytic function 
of z, with a usual branch cut for | arg z\ — n, and a singularity at ^ = 0. Therefore zo{mo,go) 
can be defined as a multivalued analytic function on the complex go plane, with a branch cut 
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for I arg go] = t: and a singularity at go = 0. So we have to consider its principal branch in the 
plane cut along the negative real axis. The analytic continuation corresponds to the definition for 
zo{mo,go) in the whole coupling constant complex plane except for a branch cut for | arg z\ = tt. 

To generate the Schwinger functions we introduce sources in the model. Thus we have that 
the zero- dimensional generating function z{mo,gQ; h) is given by 

z{mo,go; h) = sj^ dip exp (-^"^o " |y</^'') cosh{hip). (37) 

As we discussed, it is possible to find z{mo, go; h)\h=o in a closed form. Nevertheless it is not 
possible to express z{mo, go', h) in terms of known functions. If we try to expand exp{hip) in 
power series and, in order to solve the resulting integrals, we interchange the summation and the 
integration, we have problems because the power series is uniformly convergent only if \hip\ < 1. 
The result of this operation is z^^^ (mo, go', h) which is the asymptotic expansion of z{mo, go', h). We 
may write z^^\mo, go', h) ~ z{mo, go', h). Thus we have 

oo 

z«(mo,^o;/i) = E /i'Vfe(rno,5o), (38) 

fc=0 

where the coefficients are given by 



fk{mo,go)= 2A;! ^amg^^^°^"^°'^°^' ^^^^ 



Recall that zo{mo,go) is the generating function in the absence of sources. Using Eq.(36) we 
evaluate the partial derivatives of zo{mo, go) in the above formula. After some algebra, we have the 
asymptotic representation for z{mo,go', h) in terms of derivatives of the confluent hypergeometric 
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function of second kind, 

where the coefficients Ck are given by Ck — ^^/^ Let us study the singularities of 2;^^^(mo, go] h) 
in the complex coupling constant for < \go\ < oo. The derivatives of the confluent hypergeometric 
functions of second kind are given by 

^*(q;, 7; z) = + n, 7 + n; z), (41) 

where the coefficients {a)k are defined by 

(a)o = 1, ... = ^^^^^ = «(a + !)...(« + ^ - 1), (42) 

for A; = 1,2, .... Therefore, again we note that in the series representation for z^^\m,o,go; h) we 
find branch points at gQ = gQ = 00 and a branch cut at arg [go) = tt. 

An interesting question is related to the zeroes of the confluent hypergeometric function of 
second kind '^{a,c]z). Since the zero-dimensional partition function ZQ{mQ,gQ) is proportional 
to ^(a, c;^;), (a = |,c = |) it is important to know if there is a ratio between mo and g^ 
where zoirno^go) vanishes. For a and c real it is well known that '^{a,c]z) has only a finite 
number of positive zeros, and there are no zeroes for sufficiently large z [85]. Also the confluent 
hypergeometric function of second kind can not have positive zeros if a and c are real and either 
a>Oora — l-|-c>0. Consequently, there are no real values for mo and g^ where the free-energy 
per unit volume diverges. Note that \l'(a, c; z) has complex zeros for real a and c, but we are not 
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interested in such situation. In the next section we repeat our zero-dimensional analysis for the 
ultra-local quantum electrodynamics. 



5 The strong-coupling expansion and the ultra-local ap- 
proximation in quantum electrodynamics 

It is not difficult to study quantum electrodynamics also in the context of the strong-coupling 
expansion. We assume that such a theory defined in a Minkowski spacetime can be extended to 
the Euclidean formulation. Accepting this point, let us investigate the strong-coupling regime of 
the theory. The treatment is standard and the idea is the same as for the scalar theory, treating 
the kinetic terms of the theory as a perturbation and solving the self-interacting part exactly. The 
only point that we have to call the attention of the reader is that we have to introduce a mass 
term for the photon field. Here we are following the approach developed by Cooper and Kenway 
[44] and Itzykson, Parisi and Zuber [86]. 

The generating functional of all Schwinger functions for quantum electrodynamics defined in 
a d-dimensional Euclidean space is given by 



where Ai^{x) and il^{x) are respectively the gauge and fermion field and C{An,il^,il)) is given by 



Z{r],r],J^)^ / dtp dA^ exp I - / d'^x C{A^,tp,ij) + 



sources 



(43) 




(44) 
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Note that it is necessary to introduce a photon mass in the same way that we did for the case of 
the scalar theory in order to regulate the photon part of the strong-coupling expansion. Although 
the introduction of the mass term eliminates the necessity of a gauge fixing term, it is useful to 
work in a general gauge. 

The Euclidean form of the generating functional of complete Schwinger functions in a d- 
dimensional space in a covariant gauge is given by 

Z{f}, r], Jf,) ^ KaK^ J # # dAi^ exp ^- J (fx Qm^A^ + ii{im + ie-f^A^)il^ + soucesj^ (45) 

where 

= exp (i / I ,-'y^^D;ii.,y)^) , (46) 

and also 

ir, = axp(-/d^./d^,^S-'(..,)^), (47) 

where r] and fj are anticommutating Grassmann sources. 

In a covariant gauge we have D~^(a;x — y) — (S^^A — (1 — ^)dnd^^ 6'^{x — y) and also 
S~^{x,y) = 'jf_,d^S'^{x — y). As we have been discussing, the generating functional defined by 
the remaining functional integral is a product of one-dimensional integrals in each point of the 
Euclidean space. Consequently, let us study the zero-dimensional generating function of quantum 
electrodynamics defined by z{r], f],j), where we are choosing m = = 1 [87]. Thus in each point 
of the Euclidean space we have the following generating function: 

^{v,VJ) = - J dAdipdiljexp 'ip{l - eA)'ip + ipr} + fjip + jA^ , (48) 
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where ip and ip are complex conjugate c-numbers variables. Integrating over the ip and ip variables 
it is easy to show that the generating function given by Eq.(48) can be written as 

f dA ( \ 

4v,rj,j) = J (1 _ eA) "^""P i ~T ^ ^^^^ ~ ^ • ^^^^ 

Let us first study the generating functional z{'r], f], j) before using Furry's theorem and let us call 
this generating function at zero external sources as I{v), were v — —K Thus we have 

Following Stieljes [59] it is possible to show that if a function F{x) is defined by 

Fix) = r /^''\ du, (51) 
^ ' Jo (x + u) ^ ^ 

the following series J2^=i o-n " is an asymptotic expansion for the integral in which the coefficients 
a„ are given by (— l)"'~-'^a„ = /q°° f{u) vP'~^du, where n = 1, 2, ... 

This generating function for zero external sources I{z) possesses an asymptoptic expansion 

I{v) ~^a„^;-"+\ (52) 

n=l 

where the coefficients are given by (— 1)"~^ a„ = /q°^ dxx'^~^ exp(— ^x^). Using the fact that 
quantum electrodynamics is charge conjugate invariant, and making use of the Furry's theorem, 
we have that the generating function of quantum electrodynamics must be given by 



(1 - 6^2) 
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(53) 



Note that the photon propagator is given by G =< > and the electron propagator is S =< 

— > where the average is over the measure — — r exp(— 

Let us proceed in the study of the generating function for the zero-dimensional quantum elec- 
trodynamics with j, f] and rj being the c-number sources. In order to obtain some information 
about the structure of the singularities of the generating function, and also the generating func- 
tional of the Schwinger functions in the complex coupling constant plane, let first try to solve the 
integral defined by Eq.(53). It is clear that the integral that defines z{'r],f],j) is meaningful only 
for negative. To simplify our discussion let us first assume that the c-number sources fj and t] 
are zero i.e. fj = rj = 0. In this particular situation we have 

The integrand of Eq.(54) has two branch points at A = ^ and A——K To perform the integral, 
let us first make the replacement e — > ie. In this case the branch points appear in the imaginary 
axis of the A complex plane and after we impose that our function zo(ie; j) is defined only for 
> we have 

^o(ie; j) = j .,,1 exp jA ■ (55) 

•''0 {l + e^A^)2 \ ^ J 

Even after this "improvement", it is very difficult to express the integral in terms of known 
functions. Nevertheless an asymptoptic expansion for small j can be found. Using a Taylor 
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expansion for zo{ie; j) near j = we have 



zo{ie; j) = zo{ie; j)\j=o + j \j=o + ... (56) 



Let us calculate first Zo{ie; j)\j=o- We have 

Zo{ie; j)|j=o = / - — cxp(-^). (57) 
J (l + e"^/l^)2 ^ 

Thus, it is not difficult to show that zo{ie; j)\j=o is given by 

zo{ie; j)\j=o = ^ exp (^^) Ko , (58) 

where Ko{z) is the Macdonald function of zero order. It is not difficult to perform the integral 
that appears in the second term of the Taylor expansion. A simple calculation gives 



where $(x) is the Fresnel integral defined by 



$(a;) = 4= r dt exp(-t^). (60) 

VTT Jo 



Using Eq.(56), Eq.(58) and Eq.(59) we have 



Mie; S) = i (exp {±) K„ {±) + .yjexp (i,) (l - * ^ j j j + .. (61) 

It is clear that there is a branch point singularity in the complex coupling constant plane in the 
zero-dimensional model. 
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In the next section we will analyse more general scalar models with polynomial and non- 
polynomial interactions in the context of the strong-coupling expansion. Now we turn to the 
second question that we have raised: in which circumstances the strong-coupling expansion must 
be used? To clarify this problem, in the next section we perform an perturbative expansion in two 
toy-models where the use of the strong-coupling expansion is imperative. 

6 The ultra-local approximation in scalar models. 

In the infinite cut-off limit of an infrared free theory, the usual perturbative expansion where we 
assume that the non-Gaussian contribution is a perturbation of the corresponding free theory can 
not be used, and we resort to an alternative perturbative expansion. The strong-coupling expansion 
is an alternative perturbative expansion suitable for treating this situation. It is remarkable that 
in the strong coupling expansion different scalar theories can be treated in the same way, since 
we factor out the free part of the Lagrangian density and evaluate the remaining non- Gaussian 
contribution in a closed form. Prom this discussion we see that this unusual expansion can be 
performed for any polynomial or non-polynomial interaction V{gi](f), where gi,i = l,2,..n are 
the coupling constants of the model. In a different context, for the study of non-polynomial scalar 
models at finite temperature in the one-loop approximation, see for instance Ref. [88]. 

Going back, the formal representation for the generating functional of complete Schwinger 
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functions Z{h) using the strong-coupling expansion is given by 

ZW= (l-i/<i^x/d''.^ir(m„..;.-.)^ + ...) Q,(a;h), (62) 

where the ultra- local generating functional Qoio", h) is defined by the following functional integral: 

Qo((t; h)^U j [dip] exp (^j d'^x {~a ml ^\x) - V{gi; if) + h{x)ip{x))^ , (63) 

and J\f is the normalization factor. Let us study the ultra-local generating functional Qoio", h) 
in detail. Note that a naive use of a continuum limit of the lattice regularization for the ultra- 
local functional integral leads to a Gaussian theory, where we simply make use of the central limit 
theorem, which states that a probability distribution of the sum of n independent random variables 
becomes Gaussian in the limit — > oo. As ha been discussed by Klauder, in the limit where the 
coupling constant goes to zero, all of the interacting theory solutions become the pseudo-free ones. 
Then following Klauder we are able to represent Qoio", h) as 

QQ{a] h) = exp J d'^x L{a; h{x))^ , (64) 

where L{a; h{x)) is some functional. The formulae given by Eq.(63) and Eq.(64) are fundamental 
for our study. To proceed, let us see how it is possible to extract some informations selecting 
particular potentials V{gi;ip). We limit ourselves to models with only one component. The 
generalization of our investigations to models with more than one component does not present 
any difficulty. For a discussion of the strong-couphng expansion in the 0{N) model, see for 
instance [42]. 
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The ideas of the preceding sections will be ilustrated in two models where the usual perturbative 
expansion in the coupling constant can not be performed. The first one is given by interaction 
Lagrangian Cxx{P, 7; defined by 

Cxi{(5, 7; ^) = P^'ix) + 7^""(x). (65) 

where P and 7 are bare parameters and p is an integer. The second model that we would like to 
discuss, which we call the sinh-Gordon model, is defined by the following interaction Lagrangian: 

CixiiP, r,^) = P (cosh 7(^(0;) - 1) , (66) 

where /? and 7 are also bare parameters. Let us start studying the model defined by Eq.(65). Note 
that in this model we have a suppression of the configurations fluctuations around (p — 0. For p 
even the model has two minima. It is not difficult to show that the interaction Lagrangian density 
of the model has a power series representation. The equilibrium values are given by <^o = (^)^ 
and ipo = — (^)^- Let us choose the case where ipo > and define a new field (f){x) = {^{x) — (po). 

Using the binomial expansion and its generalization 

00 

(i+^r = E N<i> (67) 

A;=0 

where the coefficients of the expansion are given by 



we have 



£xx(7, 0) = E c(p, k) ct>{xf + 7 ^0 E ^(^f- (69) 

k=0 k=p+l 
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The coefficients c(p, k) are given by 

c{p, k) = [(5 il) + 7 ^o'-' t^)) ■ (70) 

Note that the generahzation of the binomial series is vahd for any complex exponent a. In other 
words, the power series in Eq.(67) is convergent everywhere in the a complex plane. Since we are 
interested in the non-perturbative effects, let us study the ultra-local version of the model. The 
zero-dimensional generating functional is given by 




It is not possible to express Z2(/3, 7; h) in terms of known functions. First let us express the 
zero-dimensional generating function Z2{/3,j; h) in the absence of sources in a closed form. Then, 
expand cosh.{hf) in power series and, in order to solve the resulting integrals, interchange the 
summation and the integration. Again 2;'^^^(/?,7; h) that we obtained after use the power series 
expansion is the asymptotic expansion of 2:2 7; h) and we write that z^'^\(5,^\ h) ~ 2;2(/3, 7; h). 
It is not difficult to find Z2{(3,^] ^)|/i=o- Using the identity 

r dx x''-^ex^^(-(5xP -^x-A = -(U^Kk {2J^\ , (72) 

JQ ^ ^ p p ^ \ ^ J 

that is valid for Re (3 > and Re^ > 0, and where K^{z) is the modified Bessel function of third 
kind, we have that the zero-dimensional generating function Z2{f3, 7; h) in the absence of sources 
is given by 



MPn; h)\h=o = 
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Using that 7; h) ~ Z2{P,'j; h) we have 



where the coefficients c(p, k) are given by 

Let us discuss now the sinh-Gordon modeL This model is also non-renormalizable in the weak- 
coupling-perturbative expansion, where f3 and 7 are the coupling constants and we choose a — 0. 
It is clear that the zero- dimensional generating function 2:3 (/3, 7; h) in the absence of sources can 
be found in a closed form and it is given by: 

zs((3,r,h)\h=o^—Ko((3). (76) 
7 

It is interesting that for this kind of theory, even in the presence of sources, we can find a closed 
form for the zero-dimensional generating fuction 2:3 (/3, 7; h), 

z,{p,r,h) = —K.{p), (77) 
7 7 

where again Kj,{z) is the modified Bessel function of third kind. Using Eq.(64), Eq.(76) and 
Eq.(77) we have that the ultra- local generating functional for the sinh-Gordon model is given by 

Qoih) =A^exp (^S'^iO) J (fx In {^Kn_{(i)^^ , (78) 

where the normalization N can be found using that (5o(/i)U=o = 1- It is remarkable that in 
the Sinh-Gordon model it is possible to extract information from a finite number of terms of the 
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infinite series representation of tlie generating functional in the strong coupling expansion. In the 
Sinh-Gordon model the ultra-local approximation is exact. Thus Z{h) = Qo{h). Let us calculate 
the "free-energy" per unit volume /(/?, 7; h) for /3 7^ and also 7 7^ 0. It is clear that we have 

fren{P,r, h) = In (^^j + In . (79) 

The first question that can be asked is related to the fact that freniP, 7; h) can diverge, since the 
modified Bessel function of third kind K^{z) has zeros in the complex plane. Prom this discussion 
one needs information about the distribution of the zeros of the modified Bessel function of third 
kind in the complex plane. For real u, K^{z) has no zeros in the region |ar5r2;| < | and in the 
complex plane with a cut along the segment (— cx), 0] the Ki^{z) has a finite number of zeros. Thus 
the free energy per unit volume in the model is finite. In the next section we sketch how it is 
possible to obtain a renormalized Schwinger functional going beyond the ultra-local approximation. 

7 Going beyond the ultra-local approximation in scalar 
models. 

In the next section we will be interested in computing global quantities, as for example the 
free energy or the pressure of the vacuum. The picture that emerges from our discussion is the 
following: in the strong-coupling perturbative expansion we reduce the problem of the singularities 
of the Schwinger functional into two parts. The first one is how to define the ultra-local generating 
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functional and the second one is to regularize and renormalize the other terms of the perturbative 
expansion, which came from the coupling between distinct points and are giving by the non- 
diagonal part of the kernel. Let us assume that a scalar field with some generic interaction 
Lagrangian is defined in a compact region of the Euclidean space. Thus it is clear that the zeta 
function regularization can be used to control the divergences of the kernel K{mo; x — y) integrated 
over the volume [89] . 

In the infinite volume limit, to complete our discussion, we have to sketch the formalism that 
can be used to obtain a regularized expression for Z^^\h) going beyond the ultra-local approx- 
imation. We are using two different regularization procedures, and it is possible to identify the 
divergent contribution in each regularized expression and a renormalization procedure is imple- 
mented with an appropriate subtraction of the singular contribution. We would like to stress that 
such kind of study in another context was performed by Svaiter and Svaiter [90]. These authors 
developed a method to unify two unrelated regularization methods frequently employed to obtain 
the renormalized zero-point energy of quantum fields. Introducing a mixed cut-off function and 
studying the analytic properties of the regularized energy as a function of the two cut-off param- 
eters it was possible to not only relate the usual cut-off method and the analytic regularization 
method, but also unify both methods. 

Let us use the ideas discussed above introducing a exponential cut-off and also an algebraic 
cut-off to regularize the kernel K{mQ, x — y). For simplicity let us assume that we have a constant 
external source, i.e., h{x) — h — constant and also u = 0. The second term of the perturbative 
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series given by Eq.(28) becomes 

Z^'\h) = I d'x I d'yKimo; x - y). (80) 

The singularities of the Z^^\h) are coming from different terms. First the second derivative with 
respect to the source of the ultra-local generating functional is singular in the continuum limt. 
Since we have been discussing how to deal with the singularities of the ultra-local generating 
functional, let us study only the divergences coming from the kernel K{mo, x — y) integrated over 
the Euclidean volume where the scalar field has been defined. In the infinite volume limit, the 
Fourier representation for the kernel K{mo] x — y) is given by 

K{mo; x-y) = -^^-^ J <fq {q^ + ml) exp {iq{x - y)) . (81) 

To evaluate the behavior of the kernel K{mo; x — y) for |x — y| small and large, as we discussed let 
us introduce two different regulators. The divergent expression given by Eq.(81) can be regularized 
using for example a exponential cut-off function /i(mo,?7; q) defined by 

/i(mo, 77; q) = exp (-77 {q^ + ml)) , Re(r7) > 0. (82) 

Another possibility is to use an analytic regularization procedure introducing an algebraic cut-off 
function f2{mo.p; q) defined by 

Mmo,p;q)^{q' + ml)'', Re(p) < - 1. (83) 

It is clear that the analytic regularization that we are using is similar to the analytic and dimen- 
sional regularization used to control divergences of the Feynman diagrams in the weak coupling 
expansion [91]. Explicit and exact integrations can be performed in both cases. 

41 



In order to carry out this program, let us study first the exponential cut-off method. The 
regularized kernel i^(mo, rj; x — y) is defined by 

1 /" 

K(mo, r];x-y) = J <fq exp {iq{x - y)) (q^ + ml) exp (-?7(g^ + mo)) . (84) 

It is clear that we can write the regularized kernel K{mo, rj; x — y) as 

1 d f 

K{mo, r]]x-y) ^ ~ (27rp ^77 J ^"'^ exp(zg(x - y)) exp + mo)) . (85) 

Since the integral in Eq.(85) is Gaussian it can be performed, and we obtain the following expres- 
sion for the regularized kernel: 

K{mo,ri; x - y) = ^J^^^-^ (v^^ exp{-riml - ^{x - yf)j . (86) 

Thus we have that the regularized kernel K{mQ,ri; x — y) can be expressed as 

K{mo,r]; x - y) ^ 

1 _ exp [-r^ml - Ux - yf] - V? + :An + 4 j • (§7) 



Note that the negative powers portion of the Laurent series expansion of K{mo, rj; x — y) around 
r] — has an infinite number of terms and the regularized expression has an essential singularity at 
a — 0. Thus, let use an alternative method, i.e. the analytic regularization procedure, that we call 
an algebraic cut-off. The same idea was presented by Kovesi-Domokos [40] in the regularization 
of the strong- coupling perturbative expansion. 
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Using the algebric cut-off function /2(mo, p; q), tlie regularized kernel K{mQ, p; x — y) is defined 

by 

K{mo, p]x-y) = J (fq exp {iq{x - y)) {q^ + ml)^+f. (88) 

The regularized kernel K{mo,p', x — y) is convergent and analytic in the complex p plane for 
Re(p) < — I — 1. As in any cut-off method we have to take the hmit p — > 0, starting from 
Re(p) < — I — 1. To perform the d-dimensional integration let us work in a d-dimensional polar 
coordinate system. Defining |a: — |/| — r and q — (?i + + •• + ?d )^ it is ^asy to show that the 
regularized kernel can be expressed in the following way: 

K{mo,p;r) = -——r- / dqq^q^ + mlY+''Jd_,{qr), (89) 
(27r)''r2~^ Jo 2 

where Ju{z) is the Bessel function of first kind of order p. Let us analyse the cases r ^ and the 
case where r — separately to make our discussion more precise: 

i) the case r = 0: this one is trivial. We have for odd d that the kernel is given by K^rno, p; x ~ 
y)\p=o — and for even d it is trivial to show that 

K{mo, p;x^y)^ _L_(^2)f (9q) 

ii) the case r 7^ 0: the more interesting case, where r 7^ can be solved evaluating the integral 
I{fj,, u; a, b) defined by 

/■oo X^~^^ 

I(li,u; a,b) ^ dx-^—^—^J^(bx), a > 0, b > 0. (91) 

t/ ix '\ a J 
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To evaluate this integral first let us start from an integral representation for the gamma function 
r{z), using the identity 

1 



1 / \ 

—— dt t" exp (-t{x^ + a^)) , (92) 



(x2 + a2)M+i r(/^ + l) 

which is vahd for Re(//) > —1 and, in order to have absolute convergence of the double integral 
that we are evaluating, we assume that the parameters /j, and u are defined in the region — 1 < 
Re(j/) < 2Re(/x) + |. Then, using an integral representation for the Macdonald's function and 

also the identity given by 

I dxx''+^exp(-a^x^) J^(6x) = -— r^-— exp -— r , a > 0, 6 > 0, (93) 
Jo {2a'^y+^ \ Aa^ J 

it is possible to show that a,b) is given by 



/(//, a,b) = 2Mr(/x+ 1) ^'^"^^"^^' 



(94) 



which is valid for —1 < Re(z/) < 2 Re(/i) + |. Using the principle of analytic continuation we have 
that the regularized kernel is given by 

K{mo,p; r) = ^ [-y) r{-l-p)-'K,^^^,{mor). (95) 

We obtained that in the limit where p — > only the r = case gives contribution to Z^^\h), since 
the Gamma function r{z) has simple poles at z — 0, —1, —2... and is an entire function of z. 
From our discussions we see that the divergence that appears in the strong-coupling perturbative 
expansion in the first-order approximation is proportional to the volume of the domain where we 
defined the fields. 
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8 The renormalized vacuum energy beyond the one-loop 
level 



Finally, in this section we would like to show that using the strong-coupling perturbative expansion 
it is possible to compute the renormalized vacuum energy of a self-interacting scalar field going 
beyong the one-loop level. Actually our results are between the one-loop and the two-loop results. 

The Casimir energy, or the renormalized vacuum energy of a free quantum field in the presence 
of macroscopic boundaries, can be derived using the formalism of the path integral quantization 
and the zeta-function regularization [89] or a different method as the Green's function method. 
In fact, these results are at one-loop level. Although, higher-order loop corrections seems now 
beyond the experimental reach, at least theoretically such corrections are of interest, as has been 
stressed by Milton [92]. 

In the present chapter we analyse the important question of the perturbative expansion of 
a self- interacting scalar field in the presence of boundaries that break translational invariance. 
Using the strong-coupling perturbative expansion we will show how it is possible to compute the 
renormalized vacuum energy of a self-interacting scalar field going beyong the one-loop level. 

Starting from the Schwinger functional it is possible to define the generating functional of 
connected correlation functions, given by lnZ{h). Consequently we have 

lnZ{h) = In (^Qo{h) - ^^^^ J d'x J d^yKinio; x - y) + ...j , (96) 

where we are assuming that the external source h{x) is constant. At the first-order approximation 
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we have 

We will concentrate only in the leading-order approximation, yielding for the generating functional 

of connected correlation functions 



\nZ{h) = --J—^^P. I d^x I d''yK{mo;x-y). (98) 



2Qo{h) dh? 

Using the results of the last section we can write the free-energy per unit volume iJi^^ as 



V 

^2, 



where f{mQ,go) is an analytic function of mo and V is the volume of the compact space. Note 
that we still have to investigate the contribution coming from the ultra-local generating functional 
Qo{h). 

Next let us use the cumulant expansion, which relates the average of the exponential to the 
exponential of the average. In our case we have the following expansion 

< >= e(<^>+5«^'>-<">')+-) (100) 

From the above equation we can see that at the first order approximation the free-energy associated 
with the self-interacting field in the presence of boundaries is given by 
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The above equation is quite interesting. We have first a contribution that is the one-loop level 
and using the heat-kernel expansion it is possible to compute the one-loop vacuum energy, but 
the ultra- local generating functional give a correction to this one-loop result. Consequently we are 
between the one-loop and the two-loop level. 

9 Conclusions 

In this paper, we first discuss the usual weak-coupling perturbative expansion and the problems 
presented by it. The weak-coupling Feynman-Dyson perturbative expansion with the respective 
Feynman diagrams is a general method to calculate the Green's functions of a renormalizable 
model in field theory. Since it is possible to express the S matrix elements in terms of the 
vacuum expectation values of products of the field operator, in particle physics it seems natural 
to perform in the majority of situations the weak-coupling perturbative expansion. Also, when 
studying for example the {X(fi'^)d model, in a four dimensional Euclidean space, since the free 
field Gaussian model gives a correct description of the critical regime when d > 4, it is natural to 
perform a perturbative expansion with respect to the anharmonic terms of the Lagrangian density. 
Nevertheless, as we discussed, there are many situations where the usual perturbative expansion 
can not be used. Consequently, we used a non-standard perturbative approach that have been 
called in the literature the strong-coupling expansion. 

From the ultra-local generating functional in different models we obtained non-perturbative 
results. We first analyse the singularities of the strong-coupling perturbative expansion in the 
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{X(p'^)d model and also quantum electrodynamics defined in a d-dimensional Euclidean space. 
Prom the discussions we can see that the divergences which occur in any scalar model in the 
strong coupling expansion fall into two distinct classes. The first class is related to the infinite 
volume-continuum hypotesis for the physical Euclidean space. The second one is related to the 
functional form of the interaction action. We showed that there is a branch point singularity 
in the complex coupling constant plane in the A<^^ model and also in quantum electrodynamics. 
Second, we discuss the ultra-local generating functional in two field theory toy-models defined by 
the following interaction Lagrangians: >Cxx(/3, 7; </?) = ^^^d the sinh-Gordon 

model {jCxii{/3,'y] ip) = /3(cosh(7<^) — 1)). A careful analysis of the analytic structure of both 
zero-dimensional partition functions in a four dimensional complex (/?, 7) space is still under 
investigation. 

Performing our expansion in bounded Euclidean volume, we sketch how it is possible to obtain 
a renormalized generating functional for all Schwinger functions, going beyond the ultra-local 
approximation. Note that the picture that emerges from the discussion is the following: in the 
strong-coupling perturbative expansion we split the problem of defining the Schwinger functional 
in two parts: the first is how to define precisely the ultra-local generating functional. Here is 
mandatory the use of a lattice appoximation to give a mathematical meaning to the non-Gaussian 
functional (actually it is not easy to recover the continuum limit). The second part is to go beyond 
the ultra- local approximation and taking into account the perturbation part. This problem can 
be controled using an analytic regularization in the continuum. Besides these technical problems 
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we still have the problem of obtaning the Schwinger functions from this approach. By summing 
over all the terms of the perturbative expansion these difficulties might be solved. We believe that 
the strong-coupling perturbative expansion is not adequate to obtain the Schwinger functions of 
the models. In the strong-coupling perturbative expansion the free energy, or any quantity that 
can be derived from the free energy, can be obtained. An interesting application is to calculate 
the free energy associated with self-interacting fields going beyond the one-loop level. 

The renormalized vacuum energy of free quantum fields has been derived using different meth- 
ods, as for example the zeta- function regularization [89]. In fact, the majority of the results in the 
literature are at the one-loop level. Although higher-order loop corrections seems now beyond the 
experimental reach, at least theoretically such corrections are of interest. In the strong-coupling 
perturbative expansion we split the problem of defining the Schwinger functional in two parts: the 
ultra-local generating functional and the perturbation part that can be controled using the heat 
kernel, zeta function regularization or any analytic regularization in the continuum. Consequently, 
a natural extension of this work in to study the renormalized zero-point energy of interacting fields 
confined in a finite volume. The strong-coupling perturbative expansion is an alternative method 
to compute the renormalized free energy or the vacuum energy associated with self-interacting 
fields, going beyond the one-loop level. This topic will be a subject of future investigation by the 
author. 
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